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I.  Introduction 


The  ion-cyclotron  instability  has  been  of  considerable  interest  in  the 
study  of  space,  as  well  as  laboratory  plasmas. Various  non-local  aspects 
due  to  magnetic  shear  were  discussed  by  Ganguli,  Falmadessc  and  the  present 
author . 5 > 6 , S , 9_  ^  novel  mechanism  for  the  generation  of  electrostatic  ion- 
cyclotron  like  modes  was  recently  proposed  by  Ganguli,  Lee  and  Palmadesso. 
According  to  this  picture,  a  non-uniform  transverse  electric  field  produces 
regions  of  negative  and  positive  wave  energies  and  a  non-local  wave  packet  can 
couple  these  domains,  with  a  resultant  transfer  of  energy  which  allows  the 
wave  mode  to  grow.  The  dispersion  relation  for  a  sharp  boundary  electric 
field  had  been  developed  which  brought  out  the  essential  physical  mechanism 
for  the  instability.  To  develop  a  treatment  for  an  electric  field  with  a 
smooth  profile  and  determine  the  ensuing  changes  in  the  dispersion  relation 
was  one  of  the  problems  to  be  explored.  Going  beyond  the  second  order 
differential  equation  formulation  was  another  generalization  to  be  considered. 
Also  to  be  considered  was  the  comparison  of  the  above  mentioned  instability 
with  the  Kelvin-Helmholtz  instability  for  the  corresponding  field  profiles. 


II.  Accomplishments 


A  brief  summary  of  the  contributions  made  during  the  grant  period  is  given 

below. 

More  realistic  continuous  velocity  profiles  Vg(x,L,a)  were  considered' 1 
instead  of  the  previous,  piecewise  continuous  (top-hat)  model'-®.  The  overall 
width  of  the  velocity  profile  is  governed  by  the  parameter  L,  while  the 
parameter  a  provides  a  continuous  transition  from  the  piecewise  continuous 
profile  (a=0),  through  a  sequence  of  increasingly  smoother  profiles,  to  a 
Gaussian.  Also  considered  were  a  trapezoid  profile  of  overall  width  L  and 
Steepness  parameter  a,  and  a  generalized  sech2(x/a)  profile.  The  Numerov 
(shooting  code)  technique  was  used  to  solve  numerically  the  second  order 
differential  equation  to  find  the  eigenmodes  for  the  Ion  Cyclotron  (IC)-like 
mode. 

Corresponding  analysis  was  also  carried  out1*  for  the  trapezoid  and 
generalized  sech^Cx/a)  profiles  to  determine  the  growth  rates  for  the 
Kelvin-Helraholtz  (KH)  instability.  An  analytical  solution  can  be  given  for  the 
trapezoid  profile. 

The  main  conclusions  were  (i)  a  magnetized  plasma  with  transverse  electric 
field  is  unstable  to  two  modes,  KH  for  long  wavelengths  (k.p £ <<  1 )  and  IC  for 
(ii)  the  KH  mode  is  important  when  ka<l,  and  dominates  for  steeper 
profiles;  (iii)  the  new  mode  (IC)  has  a  relatively  weak  dependence  on  'a'  but 
depends  strongly  on  L.  For  smooth  profiles  (sech^(x/a))  the  peak  of  the  growth 
rate  for  IC  is  comparable  to  the  KH  growth  rate,  and  the  two  instability  bands 
are  widely  separated  in  R-space.  This  mode  should  thus  be  observable  in 
simulations  or  experiments;  (iv)  since  a  sharp  profile  quickly  relaxes  to  a 
smoother  profile  (increase  in  a),  the  IC  type  mode  becomes  dominant  after  the 
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initial  transient  dynamics. 


These  results  were  presented  at  the  AGU  Meeting,  December  1986  (Ref.  11, 
copy  of  Abstract  is  attached.) 

We  have  also  considered  the  generalizations  needed  to  go  beyond  the  second 
order  differential  equation  descriptions.  Explicit  dispersion  equations  were 
obtained  for  the  4th  and  6th  order  differential  equations  level  of  description. 
Also,  ail  of  these  are  successive  approximations  to  a  full  integral  equation 
treatment.  A  full  kinetic  theory  integral  equation  treatment  was  provided  in 
Ref.  12.  Formal  reduction  leads  to  the  differential  equation  descriptions  for 
the  KH  modes  for  kp^<<l  and  the  IC-like  modes  for  kpj>l.  (Copy  of  Ref.  12  is 
attached. ) 
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Blgh  frequency  Instability  In  a  Kagnetoplasaa  for 
tarloua  Transverse  Clectrlc  Field  Profiles* 
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A  mw  Instability  Mctitrlu'  giving  rise  to  high 
frequency  short  wavelength  aeclllation*  In  •  eagneto- 


pl»sM  with  transverse  electric  field  w*s  recently 
proposed.  That  study  vu  based  on  an  Idealized 
(plece*lse  constant)  velocity  prdfild.  w#  consider 
■ore  realistic  profile*  VgU.L.a),  whe-*  variation  of 
tne  paraaeter  a  p-ovldes  a  continuous  transition  fro* 
the  piecewise  cor.  at  ant  profile  through  •  sequence  of 
Increasingly  smoother  profiles,  and  finally  to  a 
Gaussian.  The  overall  width  Is  governed  by  the 
paraaeter  L.  We  find  that  the  growth  rate  ;8 
relatively  Insensitive  to  the  variation  of  the 
steepness  pa-aeeter  a,  but  has  a  stronger  dependence 
on  the  extent  L.  Thl*  conclusion  Is  also  borne  o*t 
for  other  profile*  as  well.  W*  also  cocpare  the 
growth  rate*  for  the  fluid  7el vlr-helencl  t* 
instability  with  our  results  (considering  a  trapezoid 
profile2)  and  note  that,  as  uy  be  expected,  the 
lelvln-Helcholt*  pod*  U  Important  for  the  long 
wavelength  doaaln*.  while  the  new  Instability  la 
algnif leant  for  the  shorter  wavelengths.  1ft  fact,  for 
a  wide  rang#  of  paraaeter*.  the  two  Instability  bands 
are  quite  distinct  an  the  growth  rat*  ve-aus 
wavelength  plot.  Ta-lou*  sealings  with  physical 
parameters  are  discussed. 
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A  kinetic  theory  in  the  form  of  an  integral  equation  is  provided  to  study  the  electrostatic 
oscillations  in  a  collisionless  plasma  immersed  in  a  uniform  magnetic  field  and  a  nonuniform 
transverse  electric  field.  In  the  low  temperature  limit  (kyp,  4  1,  where  ky  is  the  wave  vector  in 
they  direction  andp,  is  the  ion  gyroradius)  the  dispersion  differential  equation  is  recovered  for 
the  transverse  Kelvin-Helmholtz  modes  for  arbitrary  values  of  kt,  where  k{{  is  the  component 
of  the  wave  vector  in  the  direction  of  the  external  magnetic  field  assumed  in  the  z  direction. 

For  higher  temperatures  ( kyp ,  >  1)  the  ion-cyclotron-like  modes  described  earlier  in  the 
literature  by  Ganguli,  Lee,  and  Palmadesso  (Phys.  Fluids  28,  761  ( 1985)  ]  are  recovered.  In 
this  article  the  integral  equation  is  reduced  to  a  second-order  differential  equation  and  a  study 
is  made  of  the  kinetic  Kelvin-Helmholtz  and  the  ion-cyclotron-like  modes  that  constitute  the 
two  branches  of  oscillation  in  a  magnetized  plasma  including  a  transverse  inhomogeneous  oc 
electric  field. 


I.  INTRODUCTION 

Shear  in  the  flow  velocity  of  a  fluid  leads  to  the  low 
frequency  and  long  wavelength  Kelvin-Helmholtz  (KH) 
instability. 1  The  velocity  shear  can  be  generated  in  a  number 
of  ways.  In  a  plasma  the  existence  of  an  inhomogeneous  elec¬ 
tric  field  component  transverse  to  the  ambient  uniform  mag¬ 
netic  field  can  provide  a  transverse  velocity  shear.  The  evolu¬ 
tion  of  the  KH  instability  in  this  configuration  has  been 
extensively  studied.2 

Recently  some  space  observations3  and  laboratory  ex¬ 
periments4  seem  to  indicate  that  ion-cyclotron-like  waves 
are  observed  for  subcritical  field  aligned  currents  and  there¬ 
fore  the  origin  of  these  waves  is  somewhat  mysterious.  A 
crucial  feature  of  these  observations  and  experiments  was 
the  presence  of  a  transverse  component  of  a  zeroth-order 
electric  field.  In  order  to  study  the  role  of  the  transverse 
electric  fields  in  the  generation  of  the  ion-cyclotron-like 
waves3,4  we  suggested  a  mechanism  based  on  the  coupling  of 
the  negative  energy  ion  Bernstein  modes  (or  the  ion-cyclo¬ 
tron  modes)  in  the  region  where  the  dc  electric  field  is  local¬ 
ized,  with  the  positive  energy  ion  Bernstein  modes  (or  the 
cyclotron  modes)  in  the  region  where  the  dc  electric  field  is 
absent.5  This  is  similar  to  the  negative  energy  wave  growth  in 
an  inhomogeneous  mirror  geometry.6  In  our  initial  theory 
we  idealized  a  typical  electric  field  profile  by  a  piecewise 
continuous  function  for  simplicity.  The  gradients  of  the  elec¬ 
tric  field  were  ignored  to  avoid  the  KH  modes  for  which  the 
second  derivative  of  the  electric  field  is  necessary.  Here  we 
use  kinetic  theory  to  obtain  the  general  dispersion  relation 
rigorously,  for  the  electrostatic  oscillations  in  a  plasma,  in 
the  form  of  an  integral  equation  for  an  arbitrary  electric  pro¬ 
file.  In  various  limits  we  reduce  the  integral  equation  to  sec¬ 
ond-order  differentia!  equations  to  obtain  the  eigenvalues. 
The  integral  equation  will  be  solved  in  a  subsequent  paper. 
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II.  THEORY 

The  equation  of  motion  of  a  charged  particle  in  a  uni¬ 
form  magnetic  field  in  the  z  direction  and  a  nortuniform  elec¬ 
tric  field  in  the  x  direction  is  given  by 

r  --= -E(x)x  +  nvxz.  (1) 

di  ~  m 

where  H  =  eBo/mc  is  the  gyrofrequency,  and  e,  m,  and  B0 
are  the  charge,  mass,  and  the  ambient  uniform  magnetic 
field,  respectively.  The  constants  of  the  motion  are  (i) 
H  —  (v\  +  u2  +  t/)/2  +  the  total  energy,  where 

E{x)  =  —  cP\f{x)/dx\  and  (ii)  Xg  =  x  +  vy/Q,  which  is 
obtained  by  integrating  the  y  component  of  (1).  Using 
vy  —  —  x)  in  the  expression  for  //we  obtain  the  Ham¬ 

iltonian  for  an  equivalent  one-dimensional  problem: 

//  =  t'2/2  +  n;(Ts  —  x)2/2  -f  (e/mWx).  (2) 

Minimizing  the  potential  of  (2)  we  obtain  the  guiding  center 
position 

£  =  *+  K  -  Ms)]/n.  (3) 

another  constant  of  the  motion  that  is  an  implicit  function  of 
Xt  and  therefore  is  not  an  independent  constant  of  motion. 
Here  V, )  =  -cEt,£)/B0. 

In  order  to  recover  the  fluid  K II  modes  as  the  fluid  limit 
of  the  kinetic  formalism,  we  will  need  to  construct  an  equi¬ 
librium  distribution  function  using  the  constants  of  motion, 
such  that  the  equilibrium  density  is  uniform  or  nearly  so. 
However,  we  would  also  like  to  be  able  to  study  the  more 
genera)  case  of  an  equilibrium  with  an  arbitrary  density  pro¬ 
file.  Therefore,  we  will  choose  the  distribution  function  to  be 
of  the  form  /0(g,f/)  =  n0(£)F0{g,H),  such  that 
fF0d}v~ const.  Then  we  obtain  an  appropriate  f0  for  the 
study  of  the  classical  KH  mode  by  setting  na( g)  =  const, 
and  an  appropriate/,  for  the  general  case  by  relaxing  this 
condition.  Such  a  distribution  function,  which  leads  to  an 
equilibrium  density  uniform  to  0(e)  for  a  constant  n0,  where 
e(  —  p,/L)  is  the  smallness  parameter,  p,  (  =  £>„/fi,)  is  the 
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ion  gyroradius,  and  L  is  the  characteristic  length  associated 
with  the  external  electric  field,  can  be  found  by  a  systematic 
procedure  and  is  given  by 

U£,H)  =  -Vexp(  -0H)g(p,  (4) 

where  jV=  nu(p  /  2jt)3/2,/?=  \/v2,v,  is  the  thermal  velocity, 
and 

g(f)  =  exp{0  [eV(£)/m  +  y2E(£)/2]}V(£)~'n, 

where  77(g)  =  1  +  y’E(£)/(l.  The  quantity  17  parametrizes 
the  magnitude  of  the  velocity  shear.  Note  that  there  are  two 
crucial  parameters  in  this  problem:  (i)  rj  and  (ii)  e.  We  will 
allow  77  to  be  arbitrary  but  positive  while  assuming  €  <  1 .  The 
equilibrium  distribution  can  be  expressed  as 

f0(£M)  =  A'  exp  [  —  (P /2)w\\ 

XexP[  -  (P/2)v]]/4W>,  (5) 

where  we  have  expanded  the  x  dependence  of  (4)  around  £ 
and  neglected  terms  of  0(e 3)  and  higher.  Here,  ui\  is 

w*  =  v\  +  vdH  +  [  V’^£)/tf]  (uy(u2y)  -  u\t  3), 

.. '  ■ ....  ;; 

where  uy  —  vy  —  (uy)  and  (  )  indicates  the  time  average.  It 
should  be  noted  that  the  simplified  zerot  harder  stateas  pro^ 
vided  in  Eq.  5,  has  a  pathological  behavior  when ‘a  becomes 
very  large.  In  our  calculations  we  use  the’ transformed  co-' 
ordinates  wl  and  <h  (see  Appendix  B)'so  that  the  simplified 
zeroth-order  state  as  given  in  (5)  agrees  with  the  exact  tc-  , 
roth-order  state  as  given  in  (4)  toordere2  but  which’  does 
not  have  this  patholopcal  pro^rty?^^1;^^^3/^^/^?^ 
Integrating  (5)  over  all  velocities  we  can  show  that  the 
equilibrium  density  distribution,  n  =  h0[l'+  ofe3)],  is  uni¬ 
form  to  order  e.  It  is  possible  to  devise  a  distribution  function 
with  density  uniform  to  any  desired  higher  order  in  e,  but 
this  is  not  necessary  here.  For  generality,  in  the  following  we 
shall  consider  a  nonuniform  equilibrium  density  profile,  i.e., 
h 0  ~  Hq(£). 

Now  using  the  definitions 

<t>(r’,t')  =exp[  —  /(<uf '  —  kyy')  ]<f>(x'),  (7a) 

<^(x')  =  j  dkx  exp (,ikxx')<fik(.k'x),  (7b) 

where  <f>  is  the  electrostatic  potential  for  the  perturbed  elec¬ 
tric  field  and  linearizing  the  Vlasov  equation,  we  obtain  the 
perturbed  distribution  function 

/lUv)  =  (J dk'x  exp[H.k'xx)]<f>k(kx) 

+  ijdk'„  <f>k (k’x)(co  -k,Vt)f  dt'A^j, 

(8) 

where  r  —  t'  —  1, 

A(r,x'yj) 

=  exp{/[& :;x' +  *,(y' -y)  +kt(z‘  -z)  -  cor]). 


rt{£)  = 


Prppfjl  d£ 


=  ^-  +  — •  (10) 
277-  77 

Here,  f„  =p/L„  and  L„  =  {[dnu{£) /d£]/nu{£)}~ '  is  the 
scale  length  associated  with  the  equilibrium  density  profile. 
In  order  to  evaluate  the  time  integral  in  (8)  we  need  the 
orbits.  The  x  component  of  ( 1 )  gives  the  net  force  in  this 
direction.  Since  vx  is  oscillatory  in  equilibrium,  (vx )  =  0. 
This  leads  to  the  expression  for  an  average  equilibrium  drift 
(uy)  =  (VE(.x))  in  they  direction.  Now  expanding  VE 
around  £  we  obtain  ' 


<»,>  =  Ve(£)  + 


Yltf) 

2flV£) 


Transforming  the  equation  of  motion  ( 1 )  into  a  frame  mov¬ 
ing  with  (vy)  in  they  direction  (i.e.,  vy  —  uy  =  vy  —  (u,)) 
and  to  the  spatial  coordinate  £  (i.e.,  x  — £)  we  obtain  the 
transformed  equations  of  motion,  ; .[ 

-JUrl-'i  , .  • 

"sjni i.i-itst.’ij’.liVhwolat  K-7»  ;’t  t  i-’i* 

tH  */  =  nV(£)uy  +  [VZ(£)/2Cl  ](<«*>  -  «J), 

tii  a ’•rjiii*  (>•>  > -  ..  7\  .  . . 

dxfri  amfi;  ait  is  ..-■■■  v-  Ol'lii-jV, 

where  we  have  neglected  terms  of  O(e^)  and  higher.  These 

:  *  i .  .  4J1  .  .  i '  •  .  t  t  1  I  f  •»  .  '  >  yr»  >;  *  *■  »  ' 

equations  lead  to  (see  Appendix  AK  * »;  <  •  v 

1  /<  >  /;»»;  .jif.-t.i.  .  •  »  1  t  .'i  .  •}  it.; 

•;r;V^w*r^s- .  .•  •  niuv 


and  u  •v.-vtew 

x'  — x  =  —  (Wi/i/^n)  [cos(>f^flr  +  d>)  —  cos4>] 

+  (t27/12772n){cos[2(d>-|-^nr)]  -cos2d>>, 

(14a) 

y'  — y  =  (iz»i/770)  [sin(<I>  +  ifrjflr)  —  sin  d>] 

—  (ix»/24775/2ri) {sin [2(0  +  t)]  —  sin  2d)} 

.  +(vy)T,  .  (14b) 

z’-z  =  vxT,  (14c) 

where 

=:  wi  +  +  yE(£){(u2y)uy  -  u]/ 3)/fl2, 

and 

u>=  yE(£)ivi/n2. 

Also  =  yfrjflt  +  4>,  where  is  the  velocity  space  angle  at 
l  =  0/  The  oscillatory  terms  of  the  order  ic  in  the  orbits  are 
not  important  except  in  the  derivation  of  the  Jacobian  of 
transformation  from  the  integration  variables  (x,t>x  ,vy )  to 
(£,U7,  ,*!>),  which  will  be  necessary  in  the  following.  For  sim¬ 
plicity,  therefore,  we  shall  ignore  the  oscillatory  terms  of 
0{w)  everywhere  except  in  the  derivation  of  the  Jacobian. 
This  restriction  can  easily  be  relaxed. 

Using  (13)  and  (14)  in  (8)  we  obtain 
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f\ U,v)  =  A/o  (£.■#)  (J  dk  '  exp(ik  'Xx)<f>k  (k'x) 

-[  dk'  <(>  (k’)\  J  ~  kyVt^^i[n'  ^  ~  a ’)  +  **£-  (^u;i/7?ft)sin«t>]} 

J  »  (w  —  Jrjn'O  —  ky{vy)  -  kxvx) 


a'  =  k[wL/Cl,  k  [2  =  k  x2/tj  +  k2y/  rj2,  a'  =  tan  '(k'x{7;.%),  g  =  x  +  uy/  fl, 
and  are  Bessel  functions.  The  projection  of  (15)  in  is 


/u,(v)=  -^~jdk'x</>k{k'x)jdxf0 


X[exp[  -i{kx  -kx)x]  -£ 


Jn-  (g)(&>  -  kyVg  ) 

^  co  —  ifrjnTl  —  ky  (vy  )  —  k.v. 


Xexp{([  (kx  -  kx  )g  +  ( n '  -  n)<S>  +  na  -  n'a']}^ 


The  perturbed  density  is  then  obtained  by  integrating  the 
perturbed  distribution  function  over  all  velocities 

nn,  =  J  dv*  dv'  du‘  J  dx  j  dk* 

X  {exp[/(A:  '  -kx)x]  -exp[/(*;  -kx)l]F}, 

(17) 

where 

F=  (co-kyV,)Y - J"l°/)JAa) - 

*'.»  co  —  n'ftjn  —  ky(uy)  —  ktv , 
Xexp[/(n'  —  7»)<t>  +  ina  —  in'a‘\.  (18) 

Equation  (17)  is  the  most  general  form  for  the  perturbed 
density  for  either  ions  or  electrons.  In  the  quasineutral  limit 
the  most  general  dispersion  relation  for  the  electrostatic 
waves,  in  the  form  of  an  integral  equation,  is  given  by 
1a  S  exp (ikxx)n°k  dkx  =  0,  where  a  represents  the  species. 

Now  we  transform  the  integration  variables  in  (17) 
from  (x,v,,vy)  to  (£,wL,< J>)  using  the  appropriate  Jacobian, 
which  in  this  case  is  (see  Appendix  B)  |/|  =  —  Jrjwp,  ex¬ 
pand  the  exponentials  in  x  and  £  around  and  retain  terms 
up  to  O(r)  to  obtain 


"n  =  - 


m  J  ~  *  J  -  •  Jo  Jo 

-  \fv 

Xexp[i(A:;  -  kx)§  )G,  (19 


-‘(kx  -kx) 


wt  cos  <t> 


-  (k'x  -kx)2 


,  W[  cos:  <1> 


We  shall  first  obtain  the  differential  equation  for  the  fluid 
KH  modes  as  the  fluid  limit  of  ( 19). 

A.  Low  temperature  but  arbitrary  shear 

Consider  =  0  and  e„  =  0.  For  a  quasineutral  plasma 
the  condition  that  the  net  perturbed  density  (19)  vanishes 
provides  the  electrostatic  dispersion  relation.  Further,  since 
A'h  =  0  the  electron  contribution  to  the  net  perturbed  density 
is  ignorahle  so  that  the  ion  perturbed  density  set  equal  to  zero 
leads  to  the  desired  general  dispersion  relation.  We  shall  now 
proceed  to  specialize  the  general  dispersion  relation  for  low 
temperature  (kypt  — 0).  Upon  <t>  integration  the  term  pro¬ 
portional  to  cos  <P  in  G  vanishes  and  n'  becomes  n  in  (18). 
Since  low  temperatures  are  of  interest  we  keep  n  =  0,  +1 
terms  in  F.  Higher  harmonics  are  associated  with  higher  or¬ 
ders  in  temperature.  The  argument  of  the  Bessel  functions 
can  be  written  as  o  =  (wL/v,)(kLL)e  and,  assuming  that 
the  factors  (kiL)  and  ( wi  /v, )  are  less  than  or  of  order  uni¬ 
ty,  we  can  expand  the  Bessel  functions  to  Off2)  so  that 


co. 


-&r- 


•») 


,  ,  2io]  k.  k ' 

2  + 


-  7/iV 


2ikyco\Cl(k  '  -  kx  )  ’ 
co)  —  7/fl2 


It  is  important  to  note  that  under  these  conditions  the  terms 
of  order  unity  in  ( 1  —  F)  drop  out.  Therefore,  the  term  pro¬ 
portional  to  ( uy )  —  Vt  —  +  p) )/2ij1~0{e2), 

which  is  responsible  for  the  KH  instability,  along  with  the 
other  terms  of  the  same  order,  are  the  leading  terms  of  (19). 
This  will  lead  to  the  dispersion  differential  equation  (22), 
describing  the  KH  modes.  However,  when  tr~0(l),  the 
Bessel  functions  can  no  longer  be  expanded  and  consequent¬ 
ly  the  role  of  the  (7(c;)  terms,  and  in  particular  the  term, 
will  diminish  in  importance.  This  situation,  which  corre- 
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sponds  10  cases  in  which  kvp,  >C>(  1 )  (f  may  still  be  small), 
may  lead  to  a  significant  change  in  the  character  of  the  mode, 
and  possibly  to  finite  gyrorac'ius  stabilization  of  the  KH  in¬ 
stability. 

In  order  to  obtain  the  dispersion  differential  equation 
for  the  KH  modes  we  need  to  obtain  nt(x)  =  /  nlkjt 
Xexp (ikxx)dkx  for  the  ions  and  set  it  equal  to  zero.  After 
carrying  out  the  wL  integration,  p]  can  be  factored  out  and 
the  resulting  equation  becomes  temperature  independent  in 
th ep,  -^0  limit.  TheJr,  integration  provides  a  delta  function, 
2n5(x  —  £),  which  makes  the  £  integration  trivial  and  con¬ 
verts  the  §  dependence  into*  dependence.  Thus,  performing 
the  s  integration  after  using  f  dk  '  <j)k  (k  '  )e\p(ik  '  £) 
=  <^(£),  we  obtain  the  second-order  differential  equation 
for  the  KH  modes  in  the  zero  temperature  limit: 


ky  V’eW 


C(*)U(x)  =0, 


where 


A(x)  =  [2kyV'E(x)(j]  +  7/'(.x)fl2] /[t>^  -  J?(x)n:], 

(23) 


B(x)  =  1  +{2[r;(x)  -  l]n2}/[w2  -i7(jc)ftJ]t  (24) 
CU)  =  n2/[(o]  -  7/(x)n2],  (25) 


andru,  =cj  —  ky  K£(x).  Equation  (22)  is  the  dispersion  dif¬ 
ferential  equation  for  the  fluid  KH  modes  valid  for  arbitrary 
shear  strength,  i.e.,  for  arbitrary  values  of  V  as  long  as 
77  >  0.  Starting  from  the  fluid  equations  one  can  also  derive 
(22).  It  should  be  noted,  however,  that  when  77  becomes 
small,  the  arguments  of  the  Bessel  functions  become  large. 
Thus  for  77  — 0  the  expansion  of  the  Bessel  functions  for  small 
argument  may  not  be  a  good  approximation  unless  the  ion 
temperatures  and  hence  p,  also  become  vanishingly  small. 
Thus  if  77  is  considered  arbitrary  then  the  fluid  theory  is  valid 
only  when  77  >  0  and  the  arguments  of  the  Bessel  functions 
are  sufficiently  small. 

Pritchett  and  Coroniti7  used  the  fluid  theory  to  obtain 


the  dispersion  equation  for  the  KH  modes  for  warm  plas¬ 
mas.  They  incorporated  the  temperature  effects  through  the 
pressure  tensor  and  arrived  at  Eq.  (16)  of  Ref.  7  as  the  dis¬ 
persion  differential  equation.  If  the  temperature  related 
terms  (which  are  included  in  the  factor  R  in  their  paper)  are 
set  equal  fo  zero  in  Eq.  ( 16)  of  their  paper,7  then  it  reduces  to 
our  Eq.  (22).  At  this  stage,  by  defining  a  transformation 
U  =  (u>,0'  —  ky{hf>)/D,  where  D  =  77  —  (aq/fl)2,  it  can  be 
shown  that  all  the  rj  related  terms  can  be  transformed  away 
and  (22)  can  be  expressed  in  the  form  of  (34),  which  repre¬ 
sents  the  classical  KH  mode  except  that  <f>(x)  is  replaced  by 
U.  Thus,  although  the  eigenfunctions  are  affected,  the  eigen¬ 
values  of  the  classical  KH  modes  remain  unaffected  by  the 
magnitude  of  the  velocity  shear  in  the  zero  temperature  lim¬ 
it.  This  is  physically  satisfying,  since  if  the  gyroradius  is  zero 
and  in  the  absence  of  an  equilibrium  x  drift,  particles  cannot 
sample  thex  direction  and  hence  cannot  experience  the  ve¬ 
locity  shear  irrespective  of  its  magnitude.  However,  for  finite 
temperatures,  the  magnitude  of  the  velocity  shear  will  play  a 
role,  as  evidenced  in  the  simulation  results7  and  the  experi¬ 
mental  results  of  Jassby.8  Unfortunately  in  the  simple  model 
for  the  temperature  assumed  in  Ref.  7  this  cannot  be  ex¬ 
plained  since  all  the  temperature  related  terms  can  also  be 
transformed  away.7  For  arbitrary  shears  the  temperature 
correction  to  the  lowest  order  involves  a  great  many  terms 
and  will  be  discussed  in  a  subsequent  paper.  Here  we  make 
the  weak  shear  approximation  but  treat  temperature  as  arbi¬ 
trary  and  find  that  the  finite  gyroradius  stabilization  of  the 
KH  modes  can  be  understood  and  predicted  by  the  kinetic 
theory. 

B.  Weak  shear  but  arbitrary  temperature 

Now  we  consider  the  limit  where  the  shear  is  not  strong 
but  temperature  and  kn  are  arbitrary.  In  this  limit,  since  the 
x  dependence  of  the  equilibrium  quantities  is  weak,  we  can 
assume  kx  ~  k  [  in  (20)  so  that  G  =  1  —  F.  Also  for  weak 
shears,  tj~  1  will  be  considered.  For  generality  we  will  also 
include  an  equilibrium  density  profile  .V(g  ).  The  perturbed 
density  after  carrying  out  the  <I>  integration  becomes 


"1  „(*> =  - 


«apa  r 


J+oc  /*  +  00  /"•  -E  oc  7*  +  x  /*  ♦  xc 

dk,  txp(ikxx)  dj  dwl  u)t  dvz  dk  x  <f>K  (k  )Na  (£)e.\p[/(A: '  -  kx  )g 

-  oc  J  —  cc  *■'0  J  —  <*  J  or. 


)  II 1  -  V 


(<y,  +  co2„ 


~qj*)JUaa)\ 
-  "tie  ~  k.v,  ) 


Here  a  denotes  the  species  and  a>  —  kyV 'ga  =<u,  +  <u,„  —  co*,  where  <o2a  =  kr  VEp*/2anda>*  =  ktea„paSla.  Inthedenomi- 
nator  we  have  approximated  0)  —  ky(Vy)  =cu,  —  <z>2,  thereby  replacing  in2  by  its  average  value  2c2,.  Now  the  uq  and  the  tq 
integrations  can  easily  be  performed  to  yield 


«>„(*)  = 


-[  dk,  exp(ikxx)  f  <f£n0(s)f  dk  *  4>k  (*  )exp[/(fc  i  -  kx  )£  ] 

,«/  -  00  J  —  «o  J  —  00 


where  T„  (b)  =  /„  ( £> )exp(  —  f>),and/„  are  the  modified  Bessel  functions  and  Z(f )  are  the  plasma  dispersion  functions.  Here 
b,  =pl(k 2  +  A- 2),  We  can  expand  T„  in  k  \p\  so  that  T„  (£, )  ~T„  (/>)  +  r '„(b)k\p2„  +  •  ••,  where  T;  =  dV„/db  and 
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b  —  k  :ypl .  Since  we  luive  consistently  retained  terms  up  to  e2,  we  neglect  the  terms  0(k\p *  )  and  higher  that  are  of  higher 
order  in  e.  Using  this  expansion  for  r„  ( b ,)  in  (27)  we  perform  the  remaining  integrals  to  obtain 


,(.x)  =  - 


[  »v  4|*u  /  v  4|*i>k.  /  dx 


where  A  2  =  mav2a/4Trn0ae2,,  is  the  Debye  length  for  the 
species  a.  The  quasineutrality  condition,  |e|(nlf  —  n,,) 
=  0,  leads  to  the  general  dispersion  differential  equation  for 
the  electrostatic  modes  in  the  low  shear  limit.  Retaining  only 
the  n  =  0  term  for  the  electrons  and  considering  p\  4p]  we 
obtain 


For  higher  temperatures  (kyp,  >  1 )  the  Bessel  functions 
can  no  longer  be  expanded  and  as  explained  earlier  the  terms 
of  the  order  e2,  such  slsco,,  become  less  important.  Neglect¬ 
ing  co 2  with  respect  toco,  and  ignoring  the  density  gradient  in 
(30)  and  (31 )  we  recover  the  starting  equations  for  the  high¬ 
er  frequency  ion-cyclotron-like  modes.5 


where 


A(x)  = 

-V  4l*,k  J 

xzH-^ — -)r"' 

V  4l*.,k  J 


Q(x)  =  1  + 


yf  CO,  +  (02  —  fel*\ 

4l*,k  J 

■  i-1—— - L  lr.(fc) 

V  41*,  k  J 

rf  1  +  (a  i  +  co2/rp  +  co*/t\ 

i  V  41*,  k  / 

r  / co i  -co2/rp\ 

A  4 1 *u  l^,  J.  ' 


Here  r  =  T,/T,,  /i  =  and  b  —  (kyp:  )2.  The  sub¬ 

scripts  on  co2  and  co*,  which  are  for  the  ions,  are  suppressed. 

In  order  to  recover  the  dispersion  differential  equation 
for  the  KH  modes  in  the  weak  shear  limit  we  set  co*  =  0, 
*,l  =  0,  and  retain  n  —  0,  ±  1  in  (30)  and  (31 )  so  that 


Further,  if  we  take  the  low  temperature  limit  [Ay/i,  <Z  I.  and 
thus  r„~  (I  -  b).  F,  ~b  /2,  r„  ■"  —  1,  I'J  -  1  j  and  the  low 
frequency  limit  (co2/fl2  4  I )  in  (32)  and  (33)  and  substitute 
these  in  (29),  we  recover  the  starting  differential  equation 
for  the  KH  modes  widely  used  in  the  literature': 


III.  NUMERICAL  RESULTS 
A.  Kinetic  Kelvln-Helmholtz  modes 

We  now  numerically  find  the  eigenvalues  of  the  general 
dispersion  differential  equation  (29),  which  is  still  in  the 
weak  shear  limit,  first  in  the  fluid  limit  to  recover  the  fluid 
results  for  the  electric  field  profile  E(x)  =  ZT0  tanh(x/Z.). 
We  then  use  the  parameters  of  the  simulation7  to  reproduce 
those  elements  of  Table  I  of  Pritchett  and  Coroniti7  that  are 
accessible  to  the  weak  shear  limit  of  the  theory.  The  differen¬ 
tial  equation  (29)  is  solved  numerically  by  a  shooting  code 
for  the  complex  eigenvalues  co.  We  assume  WKB  nature  of 
the  solutions  for  x  —  oo  and  impose  an  outgoing  energy  flux 
boundary  condition.  With  $WKB  and  <^wkb  specified  at  the 
boundaries  for  an  initial  eigenvalue  co0,  we  use  a  variable 
stepsize  integrator  to  obtain  <f>  and  <f>'  at  the  origin  where  the 
matching  condition  is  examined.  If  the  matching  condition 
is  not  satisfied  a  new  co  is  assumed  and  the  process  iterates 
until  the  eigenvalue  is  obtained  with  the  desired  accuracy. 

To  recover  the  fluid  limit  we  consider  e  =  0.1,  e„  =0, 
p  =  1837,  T’E  =  (  y°E/vt)  —  1,  u  =  *„/*,  =0.0001,  and 
r  =  3.5.  For  kyp,  =  0.02,  0.05,  and  0.08  we  obtain  yL  /  V°F 
=  0.1369,  0.1868,  and  0.1067.  These  results  coincide  with 
the  fluid  calculations  and  are  consistent  with  the  normalized 
growth  rate  against  the  dimensionless  wavenumber  plot  for 
the  fluid  KH  instability  provided  in  Fig.  2  of  Ref.  7.  Thus  the 
fluid  results  are  recovered  from  the  kinetic  dispersion  rela¬ 
tion  for  the  parameter  range  that  is  fluidlike. 

We  now  use  r  =  1 ,  to  match  the  simulation  parameters.7 
Also  since  r  =  1,  c,  =  V2u,,  where  cj  =  2(  Tr  +  T,  )/m,.  For 
VF  =v^2Ffi/c,  =  0.764  we  find  that  the  eigenvalues  of  (29) 
for  kyL  =  0.393  and  e  =  0. 19  and  0.38  are  given  by  yL  /V\ 
_=  0. 1 88  and  0. 191  while  the  corresponding  eigenvalues  for 
VF  =  1.513  are  0.184  and  0.189.  Comparing  these  with  the 
corresponding  elements  ofTable  1  in  Pritchett  and  Coroniti7 
we  find  that  our  theory  is  in  agreement  with  the  simulation 
results.7  Figures  1(a)  and  1  (b)  arc  the  eigenfunction  and  the 
profile  for  the  external  electric  field,  for  c  =  0.19  and  l'F 
=  0.764.  For  moderate  shears  (es.0.5)  better  agreement 
utth  the  simulation  can  be  expected  if  the  assumption  of 
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FIG  1  Solution  of  the  general  dispersion  differential  equation  (29)  in  the  fluid  limit.  Here  €  =  0. 19,  ktp,  —  0.074,  —  0.764,  p  =  1837,  u  =  0.0001,  and 

r  —  1.  (a)  The  real  and  imaginary  parts  of  a  typical  eigenfunction,  (b)  The  external  electric  field  profile  £(x)  =  E0  tanh(x/£)  for  e  —  0.19. 


weak  shear,  i.e.,  k  '  is  relaxed.  This  will  be  the  subject 
of  a  future  article.  Higher  shear  values,  however,  arc  inacces¬ 
sible  to  the  theory  at  the  present  simplified  differential  equa¬ 
tion  level.  For  very  high  shears  a  dispersion  relation  in  the 
form  of  an  integral  equation  as  provided  in  ( 17)  will  become 
necessary.  Also,  uncertainty  in  the  simulation  results  is  ex¬ 
pected  for  higher  shears  where  7  differs  significantly  from 
unity  (but  is  still  positive),  unless  the  initial  loading  is  in 
accordance  with  the  equilibrium  distribution  given  in  (4). 
When  shear  is  weak  and  a,  =  Vhe  [the  peak  value  of 
V’E[x)/fl]  is  small,  7  is  close  to  unity.  For  low  tempera¬ 
tures  (e  —  0).  the  equilibrium  distribution  function  (4)  can 
be  reduced  to  a  Maxwellian  shifted  in  they  velocity  by  the 
magnitude  of  the  Ex  B  drift,  which  is  approximately  VF  (x). 
Such  a  distribution  function  was  used  by  Pritchett  and  Cor- 
oniti7  for  the  initial  loading.  While  acceptable  for  small  a , 
and  especially  for  low  temperatures,  this  method  may  lead  to 
significant  relaxation  of  the  initial  state  for  large  a  „  thereby 
affecting  the  equilibrium  parameters  substantially.  Thus,  for 
large  a,,  the  interpretation  of  the  final  simulation  results7 
remains  dubious. 

We  now  study  a  different  electric  field  profile, 
E(.x)  =  Eu  seclr  (x/T).  Once  again,  to  check  the  fluid  limit 
we  use  f  —  0. 1,  r  =  1  ,/t—  1837.  and  u  =  0.0001.  Figure  2  is 
a  plot  of  the  growth  rates  and  the  real  frequency  of  the  KH 
modes  as  obtained  from  Eq.  (29)  (solid  lines)  against  b.  The 
dashed  lines  are  the  fluid  results  provided  in  Drazin  and 
Howard'1  in  their  Table  I  as  sinuous  mode.  Good  agreement 
can  be  expected  for  1.  For  larger  e,  however,  the  o)2  ap¬ 
pearing  in  the  denominator  of  (321  and  (33)  cannot  be 
treated  as  negligible.  For  larger  f,  the  denominator  co2  can 
influence  the  results  by  enhancing  the  growth  rates.  This 
important  temperature  related  effect  cannot  be  predicted  by 
the  fluid  theory  including  temperature,  as  given  by  Pritchett 
and  Coroniti.7 


A  two-dimensional  particle  simulation10  using  an  ap¬ 
propriate  initial  distribution  function  to  study  the  ion-cyclo¬ 
tron-like  modes  is  currently  in  progress.  Ultimately  we  shall 
compare  the  KH  modes  with  the  ion-cyclotron-like  modes5 
through  numerical  simulation  as  well  as  through  theory. 
Thus  we  will  us e/i  =  100,  which  is  being  used  in  our  simula¬ 
tion.10  To  study  the  finite  gyroradius  stabilization  of  the  KH 
modes  we  plot  the  growth  rates  normalized  by  the  ion-cyclo- 
tron  frequency  ft,  against  b,  for  various  values  of  u  in  Fig.  3. 
We  consider  r  =  3.5,  e  —  0.43,  and  a  mild  density  gradient 
<■,,=0.05  centered  around  —  1.33  Pi  such  that  t„ 
=  —  0.05  for  —  p,  <  x  <  +  pt  and  0  otherwise,  and 
as  the  dc  electric  field  profile  we  consider  E(x) 
=  E0  scchJ(x/L)  (see  Fig.  4).  The  growth  rates  of  the  KH 
modes  are  expected  to  reduce  because  of  the  density  gradi¬ 
ent"  but  with  the  mild  e„  used  here  this  decrease  was  not 


FIG.  2.  The  real  and  imaginary  frequencies  normalized  by  ft,  for  the  KM 
instabilities  for  the  dc  electric  field  profile  given  by  F.  =  E0  sech:  ( x/L )  are 
plotted  as  a  function  of  b.  The  solid  lines  are  the  eigenvalues  of  Eq  (29) 
while  the  dotted  lines  are  the  fluid  results  given  in  Ref  9.  Mere  (  =  0 . 1, 
r=  l,u=  1837,  and  u  =  0.0001. 
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FIG.  3.  A  plot  of  the  normalized  growth  rates  of  the  KH  modes  plotted  as  a 
function  of  6  for  a  number  of  u  values.  Here  e  —  0.43,  /i  -  100,  and  r  =  3.5. 


significant.  From  Fig.  3  we  see  that  for  a  given  e  the  growth 
rates  peak  for  a  particular  b.  The  peak  of  the  spectrum  is 
localized  for  b ~ 0.1 6  (kyp,  —0.4  and  kyL  —0.93)  and  is 
maximum  for  transverse  propagation  (u  =  0.0001 ).  As  the 
obliqueness  u  is  increased,  there  is  a  sharp  decrease  in  the 
peak  value  along  with  the  narrowing  of  the  spectrum.  Be¬ 
yond  u  =  0.0075  the  growth  of  the  KH  mode  is  substantially 
reduced. 

B.  lon-cyclotron-like  modes 

Now  we  consider  the  case  svhere  k  p,  >  1 .  As  explained 
earlier,  in  this  domain  we  can  no  longer  expand  the  Bessel 
functions,  and  consequently  the  0(e2)  terms  responsible  for 
the  KH  modes  play  only  a  minor  role.  For  the  range  in  which 
we  are  now  interested,  -here  kypt  ~3,  we  can  neglect  these 
terms  for  convenience.  Vo  explain  the  ion-cyclotron-like 
modes  we  first  resort  to  the  piecewise  continuous  field  pro¬ 
file’  (see  Fig.  5).  This  is  an  idealization  of  the  actual  field 
profile  and  we  use  it  only  to  demonstrate  the  principles  in¬ 
volved  and  to  obtain  a  good  starting  eigenvalue  for  numeri¬ 
cally  tracking  the  eigenvalues  for  a  smooth  profile.  For  a 
piecewise  continuous  profile  it  is  trivial  to  derive  the  nonlo¬ 
cal  dispersion  relation.  Setting  <o2  and  co*  equal  to  zero  in 
(30)  and  (31),  we  use  (29)  as  the  differential  equation  for 
the  modes  in  question.  In  the  region  over  which  the  electric 
field  is  localized  (we  shall  refer  to  this  region  as  region  I) 
there  is  a  Doppler  shift  in  the  frequency,  i.e., 


-  M  *•'0 


HO  4  The  equilibrium  field  and  density  configuration  used  in  the  calcula¬ 
tions  of  Fig  3  Here  —  1.33/’,.*',  -  ~-O05  when  x ^  —  p, 
x  <  *■  and  ’to  otherwise 


FIG.  5.  A  schematic  representation  of  the  piecewise  continuous  dc  electric 
field  profile. 


=  co  —  ky  VE,  while  outside  this  region  where  the 
electric  field  is  nonexistent  (region  11)  there  is  no  shift  in  the 
frequency.  This  is  the  essential  feature  distinguishing  the 
two  regions.  The  matching  condition  of  the  logarithmic  de¬ 
rivatives  of  the  solutions  of  (29)  at  the  boundary  x  —  L  /2, 
provides  the  nonlocal  dispersion  relation, 

—  *,  tan(*, /2e)  =  i/r,,,  (35) 

where  k2  =  Q^p/Aico,)  and  «•„  is  identical  to  *,  if  co,  is 
replaced  by  co.  For  details  we  refer  to  our  earlier  paper.5  The 
dispersion  relation  was  solved5  for  a  wide  range  of  param¬ 
eters  to  find  growing  modes  distinct  from  the  KH  modes.  We 
first  give  a  physical  description  for  the  origin  of  these  modes. 

The  dispersion  relation  of  the  electrostatic  ion  Bernstein 
modes  is5 


D{co)  *=  1  -  T0(b)  -  Y  r.(fr).  (36) 

n > o  co  —  n  \r 

where  k ^  —0  is  assumed.  The  energy  density  of  these  modes 
is 


U<xco 


dD 

dco 


— (l 

Vn  >0 


4^r  n2n-  n 

— r - — —  j  =  ZLi2a(£u),  cr>0. 

(co2  -rt2n-)-J 


(37) 


Clearly,  these  are  positive  energy  waves.  Introduction  of  uni¬ 
form  electric  field  in  the  x  direction  initiates  an  ExB  drift  in 
the  y  direction  and  consequently  there  is  a  Doppler  shift  in 
the  frequency,  i.e.,  co-+co,  —  co  —  kyVF.  The  energy  density 
in  the  presence  of  the  Doppler  shift  is  U'  ~coco, a  (co,),  which 
can  be  negative  provided  co,  <0.  Now  if  we  consider  the  lo¬ 
calized  field  configuration  as  shown  in  Fig.  5,  then  it  is  clear 
that  because  of  the  E  X  B  drift  the  energy  density  in  region  1 
becomes  negative  while  it  remains  positive  in  region  II.  A 
nonlocal  wave  packet  can  couple  these  two  regions  so  that  a 
flow  of  energy  from  region  I  into  region  II  will  enable  the 
wave  to  grow.  Based  on  this  simple  picture  we  can  predict 
some  gross  features  of  the  instability.  As  for  example,  using 
the  wave-kinetic  description  it  is  possible  to  obtain  the  ener¬ 
gy  balance  condition  for  the  system  from  which  important 
scalings  governing  the  growth  rate  can  be  predicted.  The 
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growth  of  the  wave  in  region  I  implies  a  loss  of  energy  from 
that  region.  By  conservation  of  energy,  t his  must  be  the  re¬ 
sult  of  convection  of  energy  into  region  II  and  any  local  ener¬ 
gy  dissipation  (5  . )  or  free  energy  release  (5  ,  )  processes  in 
region  I.  The  rate  of  growth  of  the  total  energy  deficit  in 
region  I  is  proportional  to  the  growth  rate  j*.  the  wave  energy 
density  £/,  in  region  I,  and  the  volume  of  region  I,  represent¬ 
ed  here  by  the  extent  in  the.*  direction  (Z. )  of  region  I  times  a 
unit  area  AL  in  the  plane  perpendicular  to  .*.  The  rate  of 
convection  through  A1  is  just  Vc  Uu  ,  where  Vc  is  the  group 
velocity  in  the.*  direction  and  Un  is  the  wave  energy  density 
in  region  II.  We  can  then  write  the  energy  balance  condition 
as 

y^LA^iS,  -S_  -  ycUn)Alt  (38) 

where  5+  and  S_  represent  the  source  and  the  sink  in  region 
I. The  eigenvalues  obtained  from  (35)  are  expected  tosatisfy 
the  energy  condition  (38).  For  the  situation  presently  under 
consideration  we  do  not  have  any  external  source  of  free 
energy  and  since  &n  ~0,  the  natural  dampings  are  absent. 
Therefore  S+  —  S_  =  0.  Now  it  is  clear  from  (38)  that  if  Ut 
is  negative  then  yean  be  positive  and  hence  the  growth  of  the 
wave  is  sustained  by  convection  of  energy  into  region  II  from 
region  I.  On  the  other  hand  if  U,  is  positive  then  the  convec¬ 
tion  of  energy  out  of  region  I  would  lead  to  a  negative  growth 
rate  and  therefore  to  damping  of  the  waves.  For  S  ±  =  0,  an 
important  scaling  can  be  predicted  from  (38),  i.e.,  y/Vc 
a:  \/L,  which  with  proper  normalizations  can  be  written  as 
Im (kxp, )  a.  e.  In  Fig.  6  we  plot  Iml^p, )  against  e  and  con¬ 
firm  this  scaling.  In  Fig.  6  other  parameters  arc  VE  =  2.9, 


FIG  6  A  plot  of  Im( it  ,p, )  against  f  1  he  linear  dependence  confirms  the 
scaling  y/Vn  a  \/L.  Here  l'r  =  2.9,  r  =  I,  u  =  0  0001,  and  / 1  =  1837. 


r=  l,ir=0.00Cl,//=  1837,  and  the  growth  rales  have  been 
maximized  in  b. 

We  shall  now  study  the  ion-cyclutron-like  modes  for 
smooth  profiles.  For  this  we  use  an  electric  field  profile  given 
by 

E(x)  =  E(/[A  sinh:(x/a)  +  1  ],  (39) 

where  A  =  1/sinh  2(x</a)  andx0  =  L  /  2.  Fora  —  0  (39)  rep¬ 
resents  a  “top  hat"  profile  that  reduces  to  half  of  its  peak 
value  at  x  =  x0.  As  a  increases  the  profile  becomes  smoother, 
and  ultimately,  when  a  =  a/o/sinh-1  (1),  which  makes 
A  =  1,  expression  (39)  reduces  to  E(x)  =  E0  sech2(x/o). 
The  shooting  code  used  for  the  determination  of  the  eigen¬ 
values  can  operate  best  when  the  initial  guess  for  the  eigen¬ 
value  is  not  too  far  away  from  the  actual  one.  Thus  it  be¬ 
comes  necessary  to  use  (39)  so  that  in  the  limit  a  — Owe  have 
excellent  guess  values  obtained  analytically  from  the  nonlo¬ 
cal  dispersion  relation  (35).  Starting  with  the  eigenvalue  for 
the  a  — 0  case  we  slowly  increase  a  to  track  the  eigenvalues 
for  the  profiles  with  the  desired  smoothness.  For  b~  8  we 
have  to  retain  #i  =  0,  ±  1,  ±  2,  ±  3,  ±  4,  ±  5,  ±  6  harmon¬ 
ics  and  the  associated  plasma  dispersion  functions  in  (29) 
that  are  evaluated  numerically;  thus  the  computations  for 
each  eigenvalue  is  CPU  time  intensive. 

In  Figs.  7(a)-7(d)  we  display  the  transition  of  the  elec¬ 
tric  field  profile  from  nearly  piecewise  continuous  to  smooth 
for  four  different  values  of  a.  Here  e  =  0.3.  In  Figs.  8(a)- 
8(d)  we  display  the  corresponding  wave  packets.  Other  pa¬ 
rameters  are b  —  8,  t  =  3.5,  u  =  0.01 1,/a  =  1837,  VE  —  0.6, 
jr„o  =  1.66 p„  and  e„=-0.07  if  xM  -  0.9p,  < x < 

+  0.9p,  and  zero  otherwise.  The  growth  rates  did  not  vary 
much  during  this  transition.  For  a  —  0.1,  0.707,  1.41,  and 
1 .89  the  corresponding  growth  rates  normalized  by  the  ion- 
cyclotron  frequency  y/ft,  =  0.048,  0.05,  0.037,  and  0.031. 
There  is  only  a  40%  reduction  in  the  growth  rate  from  the 
sharp  to  smooth  profile,  and  initially  in  going  from  a  +  0.1 
too  =  0.707  there  is  a  slight  increase  in  the  growth  rate.  This 
is  in  contrast  to  the  KH  instability  where  the  growth  rates 
are  dependent  on  the  second  derivative  of  the  electric  field 
and  are  therefore  very  sensitive  to  the  scale  size  variation. 

In  Fig.  9  we  provide  a  plot  of  the  growth  rate  and  the  real 
frequency  of  the  ion-cyclotron-like  modes  normalized  by  the 
ion  gy  rofrequency  as  a  function  of  b.  Here  the  profile  in  ( 39 ) 
is  used  with  a  =  1.87  and  the  rest  of  the  parameters  are  iden¬ 
tical  to  Fig.  8.  We  find  that  the  instability  is  ncaked  around 
kyp,  —  3,  which  forf  =  0.3  corresponds  to  kvL~  10  which  is 
an  order  of  magnitude  larger  than  the  corresponding  value 
of  the  peak  for  the  KH  modes.  Further  for  u  =  0.01 1  used 
here,  the  KH  modes  ere  expected  to  be  nonexistent  and  thus 
the  domain  for  dominance  for  the  two  modes  can  be  quite 
distinct.  This  contradicts  the  conclusion  in  Ref.  12  where  a 
simulation  based  on  only  one  set  of  parameters  obtained 
from  the  idealized  field  profile5  was  used  to  conclude  that  the 
KH  mode  will  always  dominate  the  ion-cyclotron-like 
modes.  Further  the  initial  loading  in  the  simulation12  (as¬ 
sumed  to  be  similar  to  that  of  Ref.  7)  is  improper  since  a,  for 
the  parameters  used  was  extremely  large  (greater  than  uni¬ 
ty),  and  consequently  the  simulation12  showed  a  strong  re- 
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laxation  of  the  initial  nonequilibrium  velocity  profile. 

In  Fig.  10  we  use  the  parameters  for  our  simulation10  ( to 
be  discussed  in  a  separate  article),  i.e.,  p  =  100,  r  =  3.5, 
e  =  0.43,  u=0.038,  e„  =  -0.05  for  xnn—p,<x<x„0 
+  Pi  and  0  otherwise  and  xM  =  1.33/a,  to  plot  the  growth 
rate  and  the  real  frequency  normalized  by  the  ion  gyrofre- 
quency.  Here  for  completeness  we  also  include  theruj  term  in 
(29)  to  compute  the  growth  rates  and  use  exactly  the  same 
dc  electric  field  and  density  profiles  as  were  used  to  produce 
Fig.  3.  The  inclusion  of  the  <y,  term  does  not  change  the 
eigenvalue  by  much.  The  peak  of  the  spectrum  is  around 
b~  14.  In  Fig.  3  we  used  the  same  parameters  to  conclude 
that  the  growth  rate  for  the  KH  modes  is  reduced  signifi¬ 


cantly  for  «>0.0075  and  the  peak  of  the  spectrum  is  around 
b~ 0.2.  Once  again  the  domains  of  dominance  for  the  KH 
and  the  ion-cyclotron-like  modes  are  quite  distinct. 

Finally  in  Fig.  1 1  we  provide  a  plot  of  the  real  and  imagi¬ 
nary  parts  of  the  eigenfrequency  co  normalized  by  the  ion 
gyrofrequency  ft,  as  a  function  of  VE,  the  peak  value  of  the 
equilibrium  Ex B  drift  velocity  normalized  by  the  ion  ther¬ 
mal  velocity.  Here  b  =  10,  r  =  3.5,  p  =  1837,  u  =  0.01 1, 
€  =  0.3,  x„o  =  1.66 Pi,  £n  =  —  0.07  when  xM-0.9p, 
<x  <xn0  +  0.9 p,  and  0  otherwise.  For  the  external  electric 
profile  we  use  (39)  with  a  =  1.87.  We  see  that  the  real  fre¬ 
quency  is  almost  linearly  proportional  to  VF,  which  is  in 
keeping  with  the  experimental  results  of  Sato  ct  at.' 
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FIG.  8.  The  real  and  imaginary  parts  of  the  corresponding  eigenfunctions  for  the  profiles  in  Fig.  7.  Here  f>  =  —  1837,  =  0.6,  r  =  3.5,  u  =  0.01 1,  and 

f,  =  -  0.07  ifx,„  —  0.9 p,  <x  +  0.9 p,  and  x^  =  1 .66/0,.  (a)  a  =  0.1,  (b)  a  =  0.707,  (c)  a  =  1.41,  and  (d)  a  =  1.89. 


IV.  DISCUSSION 

We  have  provided  a  kinetic  theory  to  study  the  electro¬ 
static  waves  that  can  be  excited  in  a  collisionless  magnetized 
warm  plasma  by  a  transverse  velocity  shear.  For  kypl  4  1  we 
recover  the  fluid  KH  modes  and  when  kypl  is  increased  we 
find  that  the  growth  rates  for  the  KH  modes  are  reduced  and 
for  large  enough  kyp,  the  KH  instability  is  completely 
damped.  Further,  the  growth  of  the  KH  modes  is  severely 
affected  by  the  parallel  dynamics.  It  seems  that  for  a  colli¬ 
sionless  plasma  the  KH  modes  can  grow  only  for  very  small 
Ar  1, .  As  kypt  becomes  of  order  unity  the  expansion  of  the 
Bessel  functions  is  no  longer  possible.  Consequently  the 
terms  of  (?(e2)  responsible  for  the  KH  modes  diminish  in 
importance.  At  this  point  the  large  kyp,  ion-cyclotron-like 
modes  dominate.  Further,  larger  kp  and  density  gradients 
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inhibit  the  classical  KH  wave  growth  while  both  these  effects 
favor  the  ion-cyclotron-like  waves. 

An  important  feature  of  the  ion-cyclotron-like  modes  is 
the  fact  that  the  real  frequency  of  these  waves  is  roughly 
around  ky  V\/l  (see  Figs.  9  and  1 1 ).  This  is  similar  to  the 
KH  waves,  and  therefore  the  two  instabilities  cannot  be  dis¬ 
tinguished  by  the  scaling  of  the  real  frequency  with  ky  VE. 

The  linear  dependence  of  the  real  frequency  of  the  ion- 
cyclotron-like  modes  on  the  dc  electric  field  was  not  explicit¬ 
ly  discussed  in  our  previous  paper.5  This  could  have  contrib¬ 
uted  to  a  misunderstanding,  which  led  Pritchett12  to 
conclude  that  since  the  modes  in  his  simulation  for  pt/L 
—  0.3  and  for  kyp,  =  0.47  and  0.94  displayed  the  linear  de¬ 
pendence  of  the  real  frequency  on  the  equilibrium  flow  ve¬ 
locity  they  could  not  be  the  ion-cyclotron-like  modes  that  we 
have  discussed.  A  similar  misunderstanding  was  also  dis- 
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FIG.  9.  A  plot  of  (he  normalized  real  and  imaginary  parts  of  the  eigenfre- 
quency  of  the  ion-cyclotron-like  instability  against  b  Here  ( =  0.3, 
u  =  0.011,  a  -  1  87,  r  =  3.5,  Vt  =0  6.  fi  =  1837.  .x„,  =  l.6^>,,  and  f, 
—  -  0  07  if -  0.9 'p,  <x>xM  +  0.9p,  and  zero  otherwise. 


played  by  Sato  et  al.A  in  discussing  their  experimental  results. 

Since  the  initial  electric  field  profile  used  in  the  simula¬ 
tion  described  in  Ref.  12  was  not  in  equilibrium,  the  system 
immediately  relaxed  [see  Appendix  A,  condition  (A21); 
hereto,  =  (p/a)  VE  =  (2.4)3  =  7.2>  1  ]  to  what  is  shown 
in  Fig.  4  of  Ref.  12.  This  is  very  different  from  the  initial 


FIG  10.  A  plot  "similar  to  Fig  9.  Here  u  —  0  038  and  the  rest  of  the  param¬ 
eters  are  identical  to  Fig  3 


FIG.  II.  The  normalized  real  and  imaginary  pa^_ts  of  the  eigenfrequency  of 
the  ion-cyclotron-like  waves  plotted  against  VE.  Here  b=  10,  r  =  3. 5, 
(i  =  1837,  u  =  0.01 1,  f  =  0.3,  a  =  1.87,  and  c„  =  —  0.07  if  —  0.9p, 
<x  <  jr„o  +  0.9 p,  and  zero  otherwise. 


profile  described  in  their12  Eq.  (2).  In  fact,  the  initial  profile 
is  characterized  by  two  scale  lengths  L  and  a  with  a  peak 
value  of  about  7>vi  while  the  final  relaxed  profile  is  more  like  a 
Gaussian  or  a  sech2(x/L)  type  characterized  by  only  one 
scale  length  L  and  with  a  peak  value  of  around  2u,  .  Also,  as 
explained, 1 2  the  spatial  extent  L  of  the  electric  field  increased 
during  the  course  of  the  simulation.  Conservatively  estimat¬ 
ing  the  broadening  of  L  by  only  20%,  and  considering  that 
the  final  profile  is  approximately  similar  to  sech2  ( x/L )  with 
p/L  =  0.25,  then  the  modes  at  kyp ,  =  0.47  and  0.94  would 
correspond  to  kyL  =  1.88  and  3.76,  respectively.  These  val¬ 
ues  will  be  larger  if  the  spatial  extent  is  broadened  more  than 
the  20%  assumed.  As  we  have  shown  through  the  analysis  of 
the  kinetic  KH  modes  and  Drazin  and  Howard’  through  the 
analysis  of  the  fluid  KH  modes  for  shear  profile  of  the  type 
sech2(x/Z,),  the  KH  modes  are  strongly  damped  for  kyL  >  1 
and  almost  nonexistent  beyond  kyL~  2.  Thus  the  mode  at 
kyPi  =  0.47  can  be  the  tail  end  of  the  KH  spectrum  but  the 
mode  at  kyp,  =  0.94  seems  to  be  completely  out  of  the  theo¬ 
retically  predicted  KH  spectrum  and  the  growth  rate  in  the 
second  mode  is  higher.  Hence,  the  conclusion  based  on  the 
simulations  of  Ref.  12  that  the  KH  mode  will  always  domi¬ 
nate  over  the  ion-cyclotron-like  modes  for  a  configuration 
with  a  localized  electric  field  perpendicular  to  an  external 
magnetic  field  is  at  best  dubious.  However,  we  do  agree  with 
the  other  conclusion12  that  the  idealized  field  profile 
(piecewise  continuous)  used  earlier'  to  demonstrate  the 
physical  principles  involved  is  not  suitable  for  simulation 
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purposes  and  that  a  strong  relaxation  from  such  a  profile  to  a 
smoother  profile  is  likely.  The  fact  that  the  piecewise  contin¬ 
uous  field  profile  is  an  idealization  was  emphasized  in  our 
earlier  paper.5  Here  we  ha\e  provided  an  equilibrium  distri¬ 
bution  (4)  that,  if  properly  loaded,  should  ensure  a  good 
equilibrium  even  for  moderate  shears.  Since  the  equilibrium 
distribution  as  provided  in  (4),  is  an  implicit  function  of  £, 
the  guiding  center  position,  it  is  not  in  a  convenient  form  for 
initial  loading  in  a  particle  simulation.  For  this  purpose  we 
will  express  (4)  in  terms  of  the  real  position  x.  From  the 
definition  (3)  for  £  we  obtain 

£-x  =  [vy  -  VE{£)]/n.  (40) 

Expanding  F£(£)  around  x  we  can  show  iteratively  that 
i  —  x  =  [v,  -  F£(x)]/?;(x) O 

-  [F"(x)/2r>’(x)n’][uy-  K£(x)]2.  (41) 

Comparing  (40)  and  (41)  we  find  that 

vf-yE(i)=  [vy-VE(x)]/r](x)+0(e2).  (42) 

By  definition  uy  =  vy  —  (uy)  and  using  (42)  along  with  the 
expression  (11)  for  (vy),  we  can  express  uy 
~  [v,  -  Ve(x)}/ii  (x)-O(e2).  Also  expanding  tj(£) 
around  x  it  can  be  shown  that  t](£)  =  ^(x)  -)-  0(e).  Using 
these  to  express  ufi  provided  in  (6),  in  terms  of  x  we  obtain 

w[  =t£  +  [vy  -  Fc(x)]V?7(x)  +0(e),  (43) 

and  therefore  the  equilibrium  distribution  expressed  in 
terms  of  x  becomes  /0  =  n0/0L  f0[ , ,  where  /0  „ 
=  -J0/2 it  exp(  —  / 3v)/2 )  and 


X(/r,(x))-,[l  -fO(r)].  (44) 

It  should  be  noted  that  when  expressed  in  terms  of  x  the  i] 
dependence  in  w]  changes  from  multiplying  the  y  compo¬ 
nent  of  the  velocity  to  dividing  it.  The  distribution  given  in 
(44)  is  the  zeroth-order  distribution  function.  The  correc¬ 
tion  to  order  e  is  given  in  Appendix  C. 

Consider  the  case  where  a,  =  V°E/Ln,  the  peak  value 
of  the  quantity  V'E/{1,  is  much  smaller  than  unity  (weak 
shear).  Now  17—  1  and  if  O(e)  corrections  are  to  be  ignored, 
then  (44)  becomes  a  Maxwellian  shifted  by  the  magnitude  of 
the  EX B  velocity  in  they direction.  Such  a  distribution  was 
used  by  Pritchett  and  Coroniti7  and  is  acceptable  for  weak 
shears  (a,  4. 1 )  especially  for  low  temperatures.  To  find  the 
correction  due  to  a,,  we  express  l/?7(x)  =  1  —  V'E/Sl  along 
with  the  assumption  that  the  temperature  of  the  system  will 
also  be  affected  so  that  /?— 0+80  such  that 
80~O(V'e/Vl).  Using  these  approximations  and 
80  =  /?K£/2flwe  can  express  (44)  as 

U  =  W /2tt )exp{  -  (0/2){v\  +  [i>,  -  F£(x)]2}) 

X(1  +  [0VE(x)/U\\{[vy  -  F£(x)]2-u2}). 

(45) 

The  correction  term  proportional  too,  was  also  discussed  in 
Ref.  7,  but  it  was  not  used  for  the  initial  loading  because  it 
was  expected  that  the  system  w'ould  make  the  necessary  ad¬ 


justments  and  that  these  would  be  small,  as  long  as  a,  -i  1. 
Thus  as  long  as  a ,  is  small  the  use  of  a  shifted  Maxwellian 
appears  to  be  acceptable,  although  (45)  describes  a  better 
initial  distribution.  For  moderate  shears,  however,  the  parti¬ 
cle  loading  should  be  in  accordance  with  (44),  otherwise 
strong  relaxation  from  the  initial  profile12  will  be  inescap¬ 
able.  Such  strong  relaxation  from  a  nonequilibrium  starting 
condition  is  invariably  accompanied  by  substantial  free  ener¬ 
gy  release,  which  leads  to  a  dynamic  state  quite  different 
from  the  quiet  equilibrium  essential  for  simulation  of  an  in¬ 
stability.  A  further  improved  initial  distribution  function 
with  the  0(e)  corrections  included  is  provided  in  Appendix 
C. 

It  should  be  remarked  that  in  most  of  the  experiments4 
and  space  observations’  there  exists  a  magnetic  field  aligned 
current  in  addition  to  the  transverse  localized  electric  fields. 
In  the  case  of  an  oblique  double  layer  the  magnetic  field 
aligned  current  can  originate  from  the  dc  electric  field  com¬ 
ponent  in  the  direction  of  the  magnetic  field  provided  there 
are  some  collisions.  As,  for  example,  in  the  experiments  of 
Alport  et  al.*  the  double  layer  has  a  component  in  the  direc¬ 
tion  of  the  external  magnetic  field  that  is  larger  than  the 
perpendicular  component,  thereby  providing  a  large  mag¬ 
netic  field  aligned  current.  Further,  in  some  recent  space 
observations”  ion-cyclotron-like  oscillations  have  been  re¬ 
ported  in  conjunction  with  simultaneous  observation  of  a 
perpendicular  component  of  a  dc  electric  field  and  a  magnet¬ 
ic  field  aligned  current  for  situations  where  the  magnitude  of 
the  magnetic  field  aligned  current  is  below  the  threshold  nec¬ 
essary  for  the  excitation  of  the  current  driven  ion-cyclotron 
instability. 14  A  recent  study 15  on  the  effect  of  the  perpendic¬ 
ular  electric  field  on  the  current  driven  ion-cyclotron  insta¬ 
bility14  indicates  that  the  perpendicular  component  of  the 
electric  field  can  lower  the  threshold  for  the  current  driven 
ion-cyclotron  instability.  The  necessary  condition  for  the 
current  driven  ion-cyclotron  instability  is  that  the  parallel 
phase  velocity  eo/k^  of  the  ion-cyclotron  waves  resonate 
with  the  parallel  electron  drift  Vd ,  such  that 

(co  —  k||  Vd)  <0.  For  subcritical  Vd,  («  —  ktl  Vd)  >0,  and 
therefore  the  Landau  damping  cannot  be  inverted.14  For 
simplicity  again  consider  the  idealized  field  profile  as  given 
in  Fig.  5.  The  introduction  of  the  perpendicular  component 
of  the  electric  field  initiates  an  ExB  drift  and  consequently 
there  is  a  Doppler  shift  in  the  frequency  co,  i.e., 
=  CO  —  ky  VE,  in  region  I  over  which  the  electric  field 
is  localized  in  the  perpendicular  direction.  Now  the  neces¬ 
sary  condition  for  the  onset  of  the  current  driven  ion-cyclo¬ 
tron  instability  in  region  I  becomes  (<u,  —  Vd )  <  0,  which 
can  be  satisfied  even  though  (co  —  kt  Vd )  remains  positive. 
Thus  the  threshold  value  for  the  magnetic  field  aligned  drift 
Vd  necessary  for  the  onset  of  the  current  driven  ion-cyclo¬ 
tron  instability  is  lowered. 

For  the  cases  where  there  is  a  magnetic  field  aligned 
drift  in  addition  to  the  transverse  localized  electric  field,  the 
term  S+  in  (38)  is  not  zero  and  can  roughly  be  estimated 
(using  the  local  theory)  to  be  proportional  to  LU, yt,  where 
the  local  growth  rate  in  the  region  1,  y,  =  —  Qx/QEa,  is 
evaluated  at  co  =  Here  QR  and  Q ,  are  the  real  and  the 
imaginary  parts  of  the  local  dispersion  relation  identical  to 
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the  expression  givei  Tn  (31)  withu>2  and  u*  set  equal  to  zero, 
and  QRu  is  the  a>  derivative  of  QR .  In  the  ion  rest  frame  the 
field  aligned  drift  Vd  provides  an  additional  Doppler  shift  in 
the  electron  term  so  that  Q ,  is  proportional  to  (ru,  —  k^  Vd  ) 
Assuming  that  the  field  aligned  current  is  also  localized 
within  the  region  I  so  that  QKal  =  L/,  /&>,  we  can  write  the 
energy  balance  condition  as 

yLA^  —  (rUj  —  )a>rLA  l  V  G  U^AL.  (46) 

We  have  neglected  the  ion  Landau  and  cyclotron  dampings. 
First  consider  the  case  where  the  electric  field  is  not  strong 
enough  to  make  e»,  <0  and  therefore  Ux  >0  but  <u,  is  less 
than  cj.  Since  cu,  <ry  it  is  possible  to  have  {io,  —  ky  VE)  <0 
when  (<y  —  ky  VB )  >  0  and  hence  the  first  term  in  the  right- 
hand  side  of  (46)  provides  a  growth  even  for  subcritical  Vd 
while  the  convection  leads  to  damping.  Now  if  <u,  <0  and 
consequently  £/,  <0,  the  convection  will  lead  to  growth  and 
the  first  term  in  the  right-hand  side  will  contribute  to  damp¬ 
ing.  However,  if  Ar(1  Vd  <0  (which  can  be  achieved  by  keep¬ 
ing  Vd  constant  but  changing  the  direction  of  parallel  propa¬ 
gation  or  vice  versa)  growth  in  region  I  may  be  expected 
from  both  the  right-hand  terms  for  (to,  —  kfl  Vd)  >0.  This 
can  be  a  likely  scenario  for  most  of  the  experiments  and  space 
observations  involving  the  ion-cyclotron-like  oscillations  for 
an  equilibrium  that  contains  a  dc  electric  field  in  addition  to 
a  uniform  magnetic  field.  More  details  will  be  provided  else¬ 
where. 


V.  CONCLUSIONS 

Using  a  kinetic  approach  we  have  studied  the  generation 
mechanisms  for  the  electrostatic  waves  in  a  magnetized 
warm  plasma  including  a  dc  electric  field  perpendicular  to 
the  external  magnetic  field.  Two  distinct  generation  mecha¬ 
nisms  are  discussed:  (i)  the  Kelvin-Helmholtz  mechanism 
and  (ii)  the  positive-negative  energy  wave  coupling  mecha¬ 
nism.  The  Kelvin-Helmholtz  mechanism,  first  discussed 
about  a  century  ago,1  depends  directly  on  the  second  deriva¬ 
tive  of  the  dc  electric  field  while  the  other  mechanism’  de¬ 
pends  on  the  inhomogeneity  in  the  energy  density  of  the 
waves.  The  KH  instability  can  dominate  for  small  kyp ,  if  the 
propagation  is  nearly  perpendicular.  For  a  collisionless  plas¬ 
ma  the  KH  instability  is  strongly  damped  even  if  k\{  is  a  tiny 
fraction  of  ky.  In  the  theory  we  have  shown  that  the  terms 
responsible  for  the  KH  wave  growth  are  proportional  to 
V  B  ( x )  and  are  of  order  e2.  Only  when  kyp,  4 1  can  the  Bes¬ 
sel  functions  be  expanded  for  small  argument  and  the  order 
unity  terms  be  dropped  out,  thereby  making  the  order  e1 
terms  primary.  This  then  assures  the  dominance  of  the  KH 
instability.  When  kyp,  is  increased  and  is  of  the  order  of  or 
greater  than  unity  the  Bessel  functions  can  no  longer  be  ex¬ 
panded  and  consequently  the  order  e~  terms  responsible  for 
the  KH  waves  cannot  gain  prominence.  At  this  stage  inho¬ 
mogeneous  energy  density  driven  modes5  dominate.  Also, 
the  dominance  of  the  KH  modes  can  be  reduced  even  for 
small  kyp,  if  more  oblique  propagation  ( larger  k(l )  is  consid¬ 
ered.  Here  we  have  also  shown  that  the  inhomogeneous  ener¬ 
gy  density  driven  modes  can  tolerate  larger  k ; .  Thus  the  two 
modes  are  quite  distinct  and.  depending  on  the  parameter 


range  (system  size,  temperature,  density  gradient,  etc.)  one 
or  the  other  can  dominate. 

It  should  be  pointed  out  that  while  the  interpretation  of 
the  inhomogeneous  energy  density  driven  modes  is  quite 
convincing  for  the  "top  hat”-Iike  profiles  as  evidenced  in 
Fig.  6,  it  is  not  so  clear  cut  for  the  smooth  profiles.  As  the 
profile  is  made  smoother,  additional  physics  is  introduced 
through  various  resonances  that  are  now  possible  since 
can  now  vary  smoothly  over  a  wide  range  of  values  as  op¬ 
posed  to  one  constant  value  in  region  I  and  a  different  con¬ 
stant  value  in  region  II  for  the  “top  hat"  profile.  Geometry 
related  effects  can  also  play  a  role.  It  was  also  noted  that  as 
the  smoothness  of  the  profile  was  increased  it  was  necessary 
to  maintain  a  very  small  amount  of  the  density  gradient  in 
the  transition  zone  in  which  the  electric  field  is  reducing  to 
zero,  to  preserve  the  growth  rates.  This,  however,  makes  the 
model  more  physcial  since  in  actual  experiments  (e.g.,  see 
Alport  et  al.4)  a  density  gradient  is  present  in  the  transition 
zone.  It  appears  that  the  density  gradient  acts  as  a  catalyst  by 
maintaining  the  growth  rate  without  much  affecting  the  real 
frequency,  although  the  exact  role  that  the  density  gradient 
plays  is  yet  to  be  fully  appreciated.  The  important  conclusion 
however,  is  the  fact  that  besides  the  KH  instability  there  is 
another  branch  that  can  also  grow  with  shorter  wavelengths 
and  higher  frequency  in  a  plasma  immersed  in  a  uniform 
magnetic  field  with  a  nonuniform  transverse  electric  field. 

Finally  we  would  like  to  point  out  that  in  the  small  kyp , 
limit  the  integral  equation  can  be  exactly  reduced  to  the  sec¬ 
ond-order  differential  equation  (22).  Thus  the  second-order 
differential  equation  level  of  description  to  study  the  nonlo¬ 
cal  wave  dispersion  properties  employed  in  this  paper  is 
more  accurate  for  the  KH  modes  than  the  ion-cyclotron-like 
modes  that  grow  for  large  kypt.  For  greater  accuracy  the 
eigenvalues  of  the  integral  equation  that  will  result  by  using 
( 17)  as  the  perturbed  density  must  be  obtained.  This  will  be 
the  topic  of  a  future  article. 
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APPENDIX  A:  PARTICLE  ORBITS  IN  A  UNIFORM 
MAGNETIC  FIELD  AND  A  TRANSVERSE  NONUNIFORM 
ELECTRIC  FIELD 

In  this  appendix  we  provide  the  derivation  of  the  parti¬ 
cle  orbits  to  Oie2 ) .  The  x  and  y  component  of  the  equation  of 
motion  ( I )  can  be  written  as 

r,=fir,-fiK{(x),  (Al) 

i,y  =  -ne„,  (A2) 

where  Ff;(.x)  =  —  cE(x)/B0.  Expressing  (Al)  in  the  guid¬ 
ing  center  framed,  and  retaining  terms  up  to  O(r)  we  ob¬ 
tain 

o,  =n(r,(£)K-  Ml)] 

-  {[IV  -  ^(£>]2/2n J}Pc(£)  +  O(r')).  (A3) 

We  now  transform  ( A2)  and  (A3)  to  a  frame  moving  with  a 
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velocity  (i\  )  in  they  direction  (i.e. 4  ((•,.))  so  that 
i\  =1  l[i)(c)ur  +  >/(i  )[<«•,■>  - 

-{[«,  +  <!',.)  -  re(£)]-’/2n-'}rj(f)).  (A4) 

An  expression  for  (r,.)  was  given  in  Eq.  (11)  in  the  text. 
Replacing  vy  by  -F  (i\  )  in  the  right-hand  side  of  ( 1 1 ),  we 
find  that  (vyj  —  TE(£)  =  TE(£)(u2)/2n77/ +  0(e~'). 
Substituting  this  in  (A4)  and  transforming  vy  to  u  +  (i>  ) 
in  (A2),  we  obtain  the  equations  of  motion  in  the  trans¬ 
formed  frames  to  0(e~): 

vx=V(£)iluy  +  [l^(f)/2n]  (<«!*> -«*),  ( A5) 

ity  =  —  Clvx.  (A6) 

Note  that  for  a  linear  field,  VE—0  and  (A5)  and  (A6) 
reduce  to  a  form  very  similar  to  that  of  the  equations  of 
motion  for  zero  electric  field  except  for  the  factor  y(£)  in 
(A5),  which  is  a  constant  for  a  linear  field.  For  V 'E  =  0  it  is 
fairly  easy  to  solve  the  equations  of  motion  and  we  can  obtain 
uy  =  A  cos  4>,  where  <I>  =  JyClT  +  <J>  and  A  is  proportional 
to  wi .  Thus  for  K£^0  we  assume  uy  —  A  cos  <t> 
-F  B  cos  20,  where  B  =  0(e*).  Differentiating  ( A6)  once 
and  using  ( A5)  for  vx  we  obtain 
«,  =  - 

=  -y(£)il7uy  -  [VE(£)/2]  (<u2>-  u7y).  (AT) 

Substituting  uy  =  A  cos  4>  +  B  cos  20,  in  the  left-  and  the 
right-hand  sides  of  (A7)  we  find 

lhs  =  —  y(£)il7A  cos  <I>  —  4y(£)fl7B  cos  2'J\  (A8) 

rhs  =  —  y(£)il 7  A  cos  <t> 

-{77(|-)n35-.[F2(^M2/4]}cos2<J>,  (A9) 
where  we  have  neglected  terms  smaller  than  0(e7).  Com¬ 
paring  the  lhs  and  the  rhs  we  find  that 

B=-VE(£)A7/  Uy(£)il7.  (A10) 

Thus 

uy  =Tcos<t>-  [  V"E  (£)A  !/12?/(£)03]cos  20,  (All) 
vx=—  —  =  4yA  sin  O - —  -  ? - sin  2<J>.  (A  12) 

n  6^n2 

The  constant  A  is  still  undetermined. 

After  multiplying  Equation  (A7)  by  iiy  it  can  be  written 
as 


w  ,  uv 


2  \  3 


(A13) 


Using  (A6)  we  can  eliminate  uy  from  (A13)  and  define  a 
constant  u>2  as 

u>l  = 

-  [F^(|)/n2][(^/3)  -wy)uy].  (A14) 

Substituting  uy  and  vx  from  (A 1 1 )  and  (A12)  into  (A14) 
and  retaining  terms  up  to  Ofr2)  it  can  be  shown  that 
A  =  w^/yjy.  Thus, 

ux  =  sin  O  —  [  FJ)(£)u;2/6j7(f),/2n2]sin  20, 

(A15) 


u y  =  -~J_  _  cos  <1>  —  [  r;  (i  )ic;/12i/(2  )2!1:  lcos  2<J>. 
xV(?) 

( A 1 6 ) 

With  the  velocities  known  it  is  a  simple  matter  to  obtain 
the  positions, 

x'  —  x  =  f  vx  dx 


=  — - — (cos  O’  —  cos  <1>) 

4y(x)0. 

V2(£)w7, 

-I - - — -(cos  20’  —  cos  2T>).  (A17) 

ny(£)2il3 

Rewriting^’  =  '  +  <f>  =  1/77  f It  -F  <J>,  where  r  =  t’  —  t 

and  <t>  =  \frjilt  -F  <t>,  we  obtain 

x'-x  =  -  [wy/JrjTgjn] 

X  [cos(yfy(£)  fir  -f-  $)  —  cos  <1>] 

+  [th/1277(^)2n] 

X(cos{2[<t>  -p  y[y(£)flT]}  —  cos  24>),  (A18) 

where  w  =  yE(£)w]/fl7.  Similarly  y'  —  y  can  also  be  ob¬ 
tained  It  should  be  noted  that  when  77 <0,  the  orbits  become 
unstable. 

For  computer  simulations  where  a  distribution  of  parti¬ 
cles  is  under  consideration  we  can  obtain  an  order  of  magni¬ 
tude  restriction  necessary  for  the  stability  of  the  orbit  of  a 
typical  particle  (characterized  by  a  velocity  v,,  the  thermal 
velocity,  and  a  displacement p,  the  gyroradius).  From  (A3) 
it  is  clear  that  as  long  as  the  first  term  in  the  right-hand  side, 
which  is  of  order  u,,  is  dominant,  the  form  of  (A3)  is 

X  =  —  n2?/(£);p  -p  corrections,  (A19) 

where  \  =  (I  —  ■*)  =  [u,,  —  VE(£)  ]/fl.  Therefore  the  re¬ 
storing  force  is  proportional  to  the  displacement.  This  en¬ 
sures  periodic  orbits  that  are  stable.  On  the  other  hand  if  the 
second  term  in  the  right-hand  side  of  (A3),  which  is  propor¬ 
tional  to  v7VE(£)/fl7,  dominates,  then  (A3)  is  of  the  form 

X  —  fl[  ^c(£)/2]t2  +  corrections.  (A20) 

The  restoring  force  is  now  proportional  to  the  square  of  the 
displacement.  Hence  the  orbits  are  no  longer  periodic  and 
therefore  become  unstable.  Thus  as  long  as  the  second  term 
of  (A3)  remains  smaller  than  the  first  term,  i.e.,  v]V'F(£)/ 
2fi2  <  y(£)v,,  we  can  expect  stable  orbits  for  the  typical  par¬ 
ticles.  This  restriction  leads  to  the  condition  for  stable  orbits 
in  a  simulation: 

p/R(£)<  [2rtv(f )/%.(£)]'«,  ( A21 ) 

where  R(£)  is  the  local  radius  of  curvature 
(  =  [  \v e  ] l/J)-  Simplifying  (A21 )  by  replacing 

the  guiding  center  position  £  by  the  real  position  x  and  con¬ 
sidering  the  electric  field  profiles  of  the  form 
VE(x)  =  V°F  f  (x/L),  wc  can  define  H(x)  such  that 

H(x)  =  [(2p,/F°)+2  e/'(x)],/J 

-e[\f"(x)/f(x)\)'n,  (All) 
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where  *  =  x/L  and  e  =  p/L.  In  order  to  have  stable  orbits 
and  avoid  (or  minimize)  relaxation  of  the  initial  electric 
field  profile  used  in  a  computer  simulation,  //(*)  should  be 
positive  for  all  x.  If  this  condition  is  violated  then  the  orbits 
will  become  unstable  and  the  profile  will  relax  until  ( A21 )  is 
satisfied. 

APPENDIX  B:  THE  JACOBIAN  OF  TRANSFORMATION 
IN  THE  VELOCITY  VARIABLES 

In  this  appendix  we  evaluate  the  Jacobian  of  transfor¬ 
mation  from  the  coordinates  (x,nx,i’y)  to  the  coordinates 
(£,wL  This  Jacobian  can  be  written  as 


dx 

dx 

dx_ 

d£ 

du\ 

<3<I> 

di\ 

dvx 

d£  »v» 

dwL 

r?<t> 

*V 

d£  t*,.-** 

dwL 

d<P 

Using  the  definition  of  £  as  given  in  Eq  (3)  in  the  text 
we  can  evaluate  the  elements  of  the  first  row  so  that 


n  da\ 

dvx 

du\ 

a", 

dwL 


n  <?<t> 

d'P 

dlL 

d<P 


where  we  have  suppressed  the  subscripts.  Multiplying  the 
last  row  by  1/11  and  adding  it  to  the  first  row  we  obtain 


v(£> 

0 

0 

dvx 

di 

du\ 

d<P 

dvy 

du„ 

dwl 

S<P 

Thus  the  determinant  has  been  considerably  simplified  and 
can  be  expanded  as 


duy 

dv,  dvy  \ 

V<9u>l 

d<P 

d'P  dwl  ) 

Recall  that  vy  =  uf  +  (vy)  and  using  the  expressions 
for  vx  and  uy  from  (A15)  and  ( A16),  and  using  equation 
(II)  from  the  text  for  {uy )  we  obtain 

vx  =  uq  sin  <I>  -  [  FE(£)a|/6r/(£),/2n2Jsin  2<I>,  (B5) 

«'i  A  ,  v'k(S)^\  f,  cos  2<J>'\ 

uy  =_— -cos<t>+  VE(£)  +  — — —  1 - —  . 

4r/(£)(l  \  3  / 


The  derivatives  necessary  in  (B4)  can  be  easily  ob¬ 
tained,  and  retaining  terms  up  to  0(e2)  we  have 

dv,  ■ 

- =  sm  '!> - — —  sin  2'I>,  ( B7 ) 

dw,  Zi/UV'-n- 


dv  wL  V'U£)w\ 

7  6 1/(^)  ft* 

dvx  Vf(£)w] 

— —  =  Wi  COS  <l> - — — -  COS  2<l>, 

dt>  377(|-)3/2n2 

avy  _  cos  $  +  Ve(£)u>l  I ^  _  cos  2<t>^ 
dwL  2j](g)2(l2  V  3  / 

Using  (B7)-(B10)  in  (B4)  we  obtain 
J  =  - 


APPENDIX  C:  EQUILIBRIUM  DISTRIBUTION  WITH 
ORDER  f  CORRECTION 

In  this  appendix  we  provide  an  equilibrium  distribution 
function  suitable  for  initial  loading  for  a  particle  simulation 
studying  the  electrostatic  waves  due  to  an  equilibrium  field 
configuration  containing  a  uniform  magnetic  field  and  a  per¬ 
pendicular  component  of  a  nonuniform  electric  field.  We 
shall  include  the  corrections  that  are  of  the  order  e  but  ignore 
the  order  e 2  corrections  for  the  time  being.  Expanding  t](£ ) 
around*  we  obtain  i7(<f)  =  rj(x)  +  (£  —  x)rj'(x)  -f  0(e2), 
where  t/' =  VE/fl.  Using  (41)  for  (£  —  *)  we  can  write 
down 

»/(£)  =  ’i(x)  +  {[vy  -  F*U)/y(.v)n]}r;;(*)/n 

+  0(e2).  (Cl) 

Also  from  the  expression  for  uy  given  in  (A16)  we  find  that 
the  time  average  (uy)  =  ui[/ 2y(f).  Expressed  in  terms  of*, 

(u2>  =  i^{l  -  [v,  -  FE(*)]F'^(*)/T,2(*)n2}/2T;(*) 

+  0(eJ). 

Using  this  and  (Cl )  in  the  expression  for  u?L  as  provided  in 
(6)  we  have 

+  Lx)]2+k  - 

V\x) 

[o,-YeW]7\  VIW 

\  2  2?7(*)  /? 72(*)fi2 


+  0{e1). 


Using  (C2)  in  the  expression  for/01  and  expanding  the  0(e) 
terms  we  obtain 

= ^ ~  f  (": + — -r*1  )i 1 1  ■  ■ 

x{i-K-^(*)][i+^(y 

7W.|+0(<a)], 

2  rj{x)  /J2j72(*)n2 


+  0(<r)  . 


IftheO(e)  term  in  (C3)  is  set  equal  to  zero  we  recover  (44). 
For  even  greater  accuracy  it  is  possible  to  obtain  the  Oie2) 
corrections  also. 
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